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Abstract 

The distribution of the Fisher zeros in the Kallen-Lehman approach 
to three-dimensional Ising model is studied. It is argued that the presence 
of a non-trivial angle (a cusp) in the distribution of zeros in the complex 
temperatures plane near the physical singularity is realized through a 
strong breaking of the 2D Ising self-duality. Remarkably, the realization 
of the cusp in the Fisher distribution ultimately leads to an improvement 
of the results of the Kallen-Lehmann ansatz. In fact, excellent agreement 
with Monte Carlo predictions both at high and at low temperatures is 
observed. Besides, agreement between both approaches is found for the 
predictions of the critical exponent ct and of the universal amplitude ratio 
A =A+/A_, within the 3.5% and 7% of the Monte Carlo predictions, 
respectively. 
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Introduction 



A proper understanding of non-perturbative phenomena in field theory and 
statistical mechanics is a major challenge in theoretical physics. Among the 
most important examples of such phenomena it counts the problem of color 
confinement in QCD, which, despite tireless attempts along the years, is still 
begging for a theoretical description. Another renowned unsolved problem is 
that of the Ising model in three dimensions. Actually, both problems are known 
to be connected, as the Svetitsky-Yaffe conjecture [2] states that the three- 
dimensional Ising model is closely related to the problem of color confinement 
and that it likely captures its main non-perturbative features near the transition. 
More generally, three-dimensional Ising model is closely related to a large class of 
physical systems near the critical point. Consequently, a deeper understanding 
of the non-perturbative dynamics of the three-dimensional Ising model would 
be of great importance in several areas of theoretical physics pQ . In turn, trying 
to find new semi-analytical methods, non-perturbative in nature, to shed new 
light on this problem represents a very interesting program. 

In a recent work [3], it has been proposed to address the problem of three- 
dimensional Ising model by using a method originally inspired in the spirit of 
Regge's theory of scattering [H [5]. Regge theory is a fully non-perturbative 
approach which allows the description of many experimental data in terms of 
an ansatz with few parameters, these to be determined in comparison with 
observations and/or by using theoretical arguments. The idea in [5J was to 
propose an ansatz, also with few parameters, to describe the free energy of 
the three-dimensional Ising model. The ansatz is heuristically motivated by 
mimicking the relation existing between the expression for the free energy in 
one and two dimensions. Then, tuning the parameters of the model for it to 
describe the high temperature regime, one ends up with an expression that also 
reproduces the results at low temperature with remarkable accuracy. 

This idea was first discussed in [9] and it was partially confirmed by the 
results of [10j . where, even in its most simplified formulation (in which a minimal 
duality breaking was assumed), such a method was shown to give results in 
a surprising agreement with observations and Monte Carlo data. In order to 
investigate the model in more detail (and with modest computational resources) 
it is necessary to find a theoretical tool able to fix (or at least to find some 
bounds on) the Regge parameters appearing in the model of [9] . Here, we report 
a remarkable progress in this direction. 

A suitable tool to theoretically constraint the ansatz of [S] is the analysis of 
the Fisher zeros, which is a powerful technique developed by M. Fisher in the 
1960's. Following the ideas of Yang and Lee [TI5], Fisher suggested to think of 
the inverse temperature j3 as a complex variable [18] . This technique permits to 
get relevant physical information about the statistical system by looking at the 
analytic properties of the (extended) thermodynamical functions. For instance, 
the distribution of zeros of the partition function in the complex plane (Fisher 
zeros) near the critical point permits to determine the universal amplitude ratio 
A = A + /A- of the specific heat; see the classical papers [2TJ [20] . It has been 
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also stressed that Fisher zeros technique turns out to be useful in analyzing 
the strength of the phase transitions (see for instance [22] [23] and references 
therein) . 

In this paper, we will be concerned with the analysis of Fisher zeros within 
the framework of [HI [TO] ■ As it will be explained in more details in the next 
sections, the ansatz proposed in [5] appears to be well suited for this kind of 
analysis since it is, in a sense, already in a Fisher zeros form. This method 
suggests a simple way to realize a strong breaking of duality, leading to an 
improvement of the semi-analytical results in comparison with the numerical 
data. 

The paper is organized as follows. In Section 1, a short introduction to 
Fisher zeros is given. In Section 2, the relations between Fisher zeros and duality 
breaking in Ising model is discussed. In Section 3, the nice interplay between 
Fisher zeros and the Kallen-Lchmann representation (KL) is analyzed in detail. 
In Section 4, we introduce what we call the strong duality breaking in the KL 
free energy for the three-dimensional Ising model. In Section 5, it is shown how 
the strong duality breaking improves the comparison with Monte Carlo results, 
both at high and at low temperatures. In Section 6, the behavior of the critical 
exponent at the critical point is described, together with the behavior of the 
universal amplitude ratio of the specific heat A. Besides, some feasible further 
improvements are pointed out. We conclude with a discussion in Section 7. 



1 Fisher zeros: a short introduction 

As it has been argued originally by M. Fisher [18], a very powerful method 
to analyze the thermodynamics of a generic spin system is to complexify its 
temperature and to study the distribution of zeros of the partition function 
in the complex j3 plane. Let Z(f3) and F(f3) be the partition function and 
the free energy of the system under analysis, respectively, and let g(r) be the 
distribution of zeros of Z(J5) in the complex /3-plane, such that g(r)dr is the 
number of zeros of Z((3) between r and r + dr, where near the critical point 
C one has /? c = (3 C + rexp(i</>) (the meaning of the parameter <\> is clarified 
below). Then, one can show that, near the critical point (3 C , the non-analytic 
part of the internal energy (that is, the singular part of the derivative of the 
free energy) and of the free energy itself, responsible for the phase transition, 
have the following form 

R 

E{(3) « dpF^f- 9{ l )dr TT-r + cc. (1) 

J /? c -/3 + rexp(z^) 
o 

R 

F w J g(r) log {(3 C — (3 + rexp(i(p)) dr 
o 

where "c.c." means the complex conjugate and R is a suitable cut-off. The 
angle <f> is the angle formed by the tangent to the curve of zeros of Z(j3) and 
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the real axis at the critical point (3 C . From the formulas above, one can derive 
an expression for the singular part of the specific heatQ 




C = -Pt \ I .. 3{r)dT . „ 2 +c.c. 







(/? c -/3 + rexp(i0)) 



Then, let us consider a system undergoing a second order phase transition, 
so that the singular part of the specific heat reads 

C = A + \P-P c \~ a 0<0 C , a>0, 

C = A_\p-p c \~ a 0>f3 c ; 

C = A + \og\0-0 c \ 0</3 c , a = 0, 

C = A_ log |/3 - /3 C | 0>P e , 

where a is the critical exponent and, in general, A + ^ A_ . It is worth mention- 
ing that the ratio A = A + /A_ is a universal quantity, so that its computation 
is a relevant (and challenging) question. The ratio A can be measured in many 
interesting physical systems. It can be shown (see, for instance, [3T] [5D]) that 
both g(r) near the critical point and the value of A are related to <f> in such a 
way that, for a ^ 0, the relations read 

a _ cos[(2 - a)4>] !_ Q 



cos[(2 - a)4> + i 



g(r) fa g r 



so that the amplitude ratio is trivial (namely A = 1) when the intersection 
between the zeros and the real axis is vertical (corresponding to <j> = 7r/2). In 
the exact solution of the Ising model in two dimensions 25j A = 1, but it is not 
the case in three dimensions, where it is known that 

A 3D » 0.55 , cb 3D ? vr/2. 

(see, for an updated review, pQ.) 

The aim of the present paper is to discuss the relation between the existence 
of a cusp in the Fisher curve (i.e. the fact that A / 1) and a breaking of the 
2D Ising self-duality. Such a duality symmetry is absent in the 3D Ising Model; 
consequently, any suitable attempt to dimensionally extend analytic methods 
from two to three dimensions would require to specify a precise mechanism of 
self-duality breaking. In [10] self-duality was broken in the softest possible way. 
Here we will show that a stronger breaking is favoured by the comparison with 
observations. 



1 here we will only consider the case of second order phase transitions 
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2 Fisher zeros and duality breaking 



Let us shortly describe the case of the two-dimensional ferromagnetic Ising 
model. The free energy is given by 



7T 

F 2D (I3) = log 2 cosh 2/3 + ±- J dtlog j 1 1 + ^1 - fc 2D (/3) 2 sin 2 t j 



(fc 2 D(/3))' 



cosh 2/3 coth 2/3 



, exp(2/3) - exp(-2/3) 
: (exp(2/3)+exp(-2/3)) 2 



(2) 



(3) 



< (fc 2 D(/3)f<l. 



The critical point of the two-dimensional Ising model is located at the max- 
imum of &2.d(/3), namely 

{k 2D {P c )f = 1 
and, being k 2 ]j & smooth function, it is also true that 



= 0. 



The function k 2 D (a sort of effective coupling constant) encodes the duality 
properties of the model: if one writes k 2 o in terms of r = tanh (3 (the high 
temperatures variable) or in terms of u = exp(— 2/3) (the variable at low tem- 
peratures) it turns out that k 2 D looks the same, and this is a convenient way to 
express the self-duality of the 2D Ising model (the self-duality of the 2D Ising 
model was discovered in [23] , before the discovery of its exact solution in 1944). 
From the exact solution (HI)-© one can determine the Fisher zeros of the 2D 
Ising model, which lay on a circle in the complex u = exp(— 2/3) plane. The 
equation determining the curve on which the Fisher zeros lay is: 



(k 2D (/3)Y 



= 1 



where /3 is now a complex variable; while in the the complex u-plane Fisher 
zeros have the following form: 

u = 1 — \/2 exp(i9) . 

Being the Fisher zeros located on a circle whose center is on the real axis, the 
intersection of the Fisher zeros with the real axis is vertical and thus the angle </> 
is it/2, implying A = 1. In the two-dimensional Ising model, the reason why the 
intersection at the critical point is vertical is simple: the physical singularity 
of the two-dimensional Ising model is the real positive solution (in /3) of the 
equation 

l-(fc 2 Z3(/3)) 2 (sin;) 2 = (4) 
which cannot be fulfilled unless 



(sin*) 2 = (k 2D (f3)f 
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Then, let us call the complex zeros near the critical point (3* — (3 C + x + iy 
(near the critical point x and y are both small). Thus, x is implicitly defined in 
terms of y by the equation 

H 2D {x,y) = 1 - \k 2D (j3 c + x + iyj 2 \ = 0. (5) 

Near the critical point, when the intersection is not vertical, one can think 
of X clS cL monotonic function of y (or viceversa). Therefore, in order for x to 
be defined implicitly in terms of y by the equation above near x = y = 0, the 
derivative of H2D{x,y) with respect to y should be different from zero. On the 
other hand, in order to have only one physical singularity, one has to demand 
the existence of only one solution of Eq. ^ on the real positive /3-axis and such 
a solution necessarily occurs at the maximum of k\ D . 

As long as k\ D (J3) is a smooth function that admits an analytic extension 
near (3 C , the derivative of H 2 D(x,y) (defined as in Eq. ([5])) with respect to 
y vanishes at x = y = 0, and this represents an obstruction to think of x as 
implicitly defined in terms of y through ([5]) at that point. So, no cusp in the 
Fisher curve could exist since the curve of distributions of zeros intersects the 
x-axis vertically. 

On the other hand, in physics one is often interested in systems with non- 
trivial ratio A / 1, where the cusp manifests itself. For instance, this is what 
happens in the Ising model in three dimensions, in what we are interested here. 
In order to analyze the possibility of achieving a realization yielding A ^ 1, first 
we have to make some comments on the general form of the partition function: 
let us assume for a moment that for a systems with only one phase transition, 
near the critical point, the equation that determines the singularities can be 
written in the form 

l-(fc e //(/?)) 2 (sini) 2 = (6) 

for some real function k e ff{(3) whose only positive maximum occurs at the 
critical point /? c , being t some dummy integration variable as in (JTJ) . In such 
a case, there is only one positive real solution occurring when both (fc e //(/?)) 
and (sin£) 2 are at their maxima (equal to one). 

Thus, Fisher zeros can be described by an equation which is formally anal- 
ogous to the one of the two-dimensional Ising model; namely 



1 



(keff(P)) 2 ■ (7) 



This is part of the proposal in [9] . Actually one may introduce at least implicitly 
a suitable fc e / / for many statistical system (with standard Fisher zeros) in such 
a way that ([7]) represents precisely the Fisher zeros of the system of interest. 

Equations © and ([7]) simply encode the fact that there is only one physical 
singularity. Function k e f f (/?) indicates how far from being self-dual a system is. 
Therefore, according the discussion above, having a value A ^ 1 would imply 
k e ff(f3) in ([6]) to be non-differentiable at the critical point, where it takes its 
maximum value. Actually, in order to achieve a cusp in the distribution of the 
Fisher zeros, it is enough to require a non-analiticity of k e ff(f3) at the critical 
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point. As a matter of fact, our results indicate that the first derivative of fe e / / (/3) 
is continuous while the second derivative is discontinuous. This is one of the 
hints for constructing our ansatz, as we are reminded of the fact that in three 
dimensions A 3 d < 1. 

In [TO], for the sake of simplicity (and because of the computational re- 
sources) a minimal duality breaking was assumed, according to what k e ff(0) 
still was smooth at the critical point. As a result, the agreement with Monte 
Carlo data and observations was found to be very good. Nevertheless, as it will 
be explained in the next sections, when the cusp (A ^ 1) in the Fisher curve is 
implemented, the comparison with Monte Carlo data turns out to be even bet- 
ter, showing a new and more efficient way of exploring the space of parameters 
of the model. 



3 Fisher zeros in the Kallen-Lehmann represen- 
tation 

A powerful non-perturbative technique in field theory is the Kallen-Lehmann 
representation. Such a tool allows to encode in a very natural way many analyt- 
ical properties of the spectral functions. It is interesting to adapt this method 
(often exclusively associated to Quantum Field Theory) to the analysis of the 
3D Ising model. This representation [3] gives rise to an ansatz for the free energy 
of the three-dimensional Ising model of the following forrr0 



F& X) (f3) = F 2D (P) + 



log 



A 



7T 7T 



(27T) 



dz dy ■ 



o o 















I* 







, (A(z)-l) Cl 
' A(z) 



sin 2 y 



where 



A(z) = (l+(l-fc e// (/3) 2 sin 2 z)< ) 2 , Co>Ci,C 2 ,C 3 >0, 



< (k eff (P)) 2 <1, l<A(z)<4, 



(8) 



(9) 
(10) 



and where the values of the parameters Coj Cii C2J C3 (henceforth Regge Parame- 
ters) in the case of the two-dimensional Ising model would correspond to A = 1 , 
Q = 1/2 fori = 0,.., 3. 

Expression ijHJ) will be referred as Kallen-Lehmann free energy. The heuristic 
argument that suggests to write down © goes as follows: think of an operator 
Oi^2 which "dresses" the trivial one-dimensional solution of the Ising model 
giving rise to the Onsager solution. Then, the Kallen-Lehmann ansatz for the 



2 The "Regge parameters" appearing in these formulas are related with those appearing in 
[9] by the following identities ( ± = u, f 3 = 1/2 = f 2 > Co = a - 



7 



free-energy for the three-dimensional Ising model is obtained by modifying Oi^2 
in such a way that the parameter^ A = 1 and ( i = 1/2 now become free 
parameters, and then applying such a modified dressing operator 02^3 to the 
Onsager solution. The result is Eq .([5]). 

A nice feature of this dimensional-recursive method is that it permits to 
distinguish between the problem of fixing five external Regge parameter^ A 
and (i = 0, ..,3), and the problem of determining the function fc e //(/3). The 
function fc e //(/3) plays the same role that £; 2 d(/3) plays in ©. Namely, it encodes 
the information about how different the degrees of freedom are in the low and 
the high temperature regimes. Thus, at least in principle, one could fix k e ff((3) 
by analyzing how the duality is broken in the three-dimensional Ising model. 
Then, Regge parameters can be found by looking, for instance, at the high 
temperatures behavior. 

An important piece of information is given by knowing which are the more 
relevant parameters in writing down the ansatz ([8|), as well as understanding 
their physical interpretation. In particular, as it will be discussed in the next 
sections, one learns that tuning appropriately just three Regge parameters is 
sufficient to get an excellent agreement, both at high and at low temperatures. 
Besides, a very good description near the critical point is also observed. Our 
results suggest that the agreement with observations and Monte Carlo data 
improves its performance by resorting to a careful analysis of duality breaking, 
instead of trying to select the Regge parameters ad hoc. From our analysis 
herein we achieve a substantial improvement of the method of [9]. 

Nevertheless, it is worth mentioning that our analysis is still far from being 
an exhaustive exploration of the whole space of parameters of the model. In 
fact, this is what makes the results more surprising, as the agreement with 
Monte Carlo data is excellent despite only a small piece of the moduli space was 
analyzed. 

Let us now compute the internal energy in the Kallen-Lehman representa- 
tion. A simple computation yields 



3 Namely, in the two-dimensional Ising model, the parameters which appear in Oi— ,2 are 
A = 1, = 1/2. In order to change the critical exponents one needs to change (at least) A 
and f 4 . 

4 In |10p . for sake of simplicity, the parameters and f 3 were fixed to 1/2; that is, to 
the values they take in the two-dimensional problem. Thus, already by using only three 
parameters, very good results in comparisons with Monte-Carlo data were obtained. 



e: 













(ii) 
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where 



N = 1 



, (A(z)-l ) 
' A(z) 



sin 2 ?/ 



N > 1. 



D = l-R c ? 



sin 2 y, 



i? 



,(A(z) 



A(z) 



(12) 



and 



<9A 



a(fce//(/3) s 



(A(z) 



-2(l+(l-fc e //(/3) 

{[2+(l-*e//C8) 



2-2 \ Cc 

sin z 



(Co 



sin 2 z) 



(l-keff([3) 2 



2 sin 2 z) Co 



(/3) 



2 sin 2 z" Co 



)1} 



<9 A i? = 2 



(A(z) - 1)^ 



[l + td-l^z)] 



It is worth noticing that, from the point of view of Fisher zeros distribution, 
there is a special value for In order to have a "singularity equation" similar 
to the one arising in the two-dimensional Ising model, it is necessary to ask that 
all the singular terms of the internal energy have their origin in factors like 

(l-fc e// (/3) 2 sin 2 zf (13) 

where the exponents ei depend on the Regge parameters £j. A potentially 
disturbing term in (fTT|) . which is not of the form (fT3|) . is iX 3-1 , since D = 
sin 2 y could be negative and (3 — 1 could be negativ^]. Then, to avoid 
undesired divergences in the internal energy, one may ask whether a special of 
Ci exists such that the maximum^ too of R 

,(A-ir' 



mn = max R = max 

A A 



A 



is less than or equal to one, and such that Too does not depend on z and f3: 

,(A-1) C1 - 



< 



< too = 1. 



Interestingly enough, such a special value of ?i does exist and it is precisely 
1/2, the value of the two-dimensional Ising model. Thus, with q\ = 1/2, when 
function R in (fT2"|) attains its maximum too (this happens for A = 2), it is always 
multiplied by its derivative 8aR (see Eq. (fT3"]) ). which vanishes. In this way, all 
the potential singular terms in Eq. (jlip are of the Fisher form in Eq. (fl"3"|) so 
that they can be analyzed with the method described in the previous sections. 
The condition £1 = 1/2 (which in [lOj was imposed for simplicity) is obtained 
here as being the appropriately value from the point of view of the Fisher zeros. 
Thus, the Kallen-Lehmann free energy is a smooth function everywhere apart 
from the critical point singled out by the equation ([7]) . 

5 This term is also disturbing since it could give rise to complex numbers in the case D is 
negative and £3 — 1 is negative. 

6 Where one "maximizes" C over the A fulfilling Eq. (1 1 1 . 
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4 Duality breaking and k e ff{{3) 



The previous analysis, relating Fisher zeros to duality breaking, suggests how 
to break duality taking into account the non-triviality of the amplitude ratio 
A 7^ 1 in three dimensions. A second important ingredient in the discussion is 
the Marchesini-Shrock symmetry [26 , which states that the partition function 
of Ising model on regular hypercubic lattices is invariant under 



(3- 



P + in—, n € Z. 



The simplest possible choice of fc e //(/9) which satisfies this constraint is 

d 3 exp(2/3) - d 2 exp(-2/3) 



keff(P) 



(di exp(2/3) + d exp(-2/3)) z 



In [10], two of the parameters di (i — 0, .., 3) have been fixed by asking the 
expected transition to occur at 

k eff (/3 c ) = l, /3 C = 0.22165, (14) 

and by requiring the vanishing of the derivative of k e ff(f3) in (3 C , namely 

P*\d k eff (P)\ = *=O (15) 

On the other hand, the other two parameters di were chosen to be equal 
to those values they take in the two-dimensional model. However, from the 
previous discussion, we learn that condition (fT5|) should not be kept in the cases 
of systems exhibiting non-trivial amplitude ratio A. Thus, we will assume that 
j3 c is a maximum of fc e // but with a left-derivative different from the right- 
derivative at f3 c . 

Indeed, there are an infinite number of possible choices of parameters that 
leads to strong duality breaking such that fc e //(/3) has a cusp at the critical 
point. As in [TU], we will follow a simplicity criterion: we will assume that both 
for (3 < [3 C and for (3 > (3 C the function fc e //(/3) has a form similar to that of 

{(3) in two dimensions (the numerical results indicate that the discontinuity 
appears in the second derivative, while the first derivative at the critical point is 
continuous). In addition, we will consider two constraints on the parameters: the 
first is that, for (3 < f3 c , fc e //(/3) is decreasing, while for (3 < (3 C it is increasing; 
the second is the continuity of fc e //(/3) together with Eq. (fT4]) . Summarizing, 
we have 



kef f (13) 

keff(P) 



4 +) exp(2/3) - 4 +) exp(-2/3) 



{d[ +) exp(2/3) + 4 +) exp(-2/3) 
4 _) exp(2/3) - 4 _) C M~ 2 P) 



4 _) exp(2/3) + 4 5 cxp(-2/3) 



,/3 < /? c (16) 



,/3>/3 c (17) 



10 



obeying 

= fc$(0c) = l> (18) 

and 

d k ef f(P) > , 0<0 e , (19) 

^fc e //(/3) < , 0>0 B . (20) 

Eq. (US)) fixes one of the d[ + ^ (as well as one of the d^) in terms of the 
others and the critical temperature, while Eqs. (fT§|) and (f2"0"|) ensure that /3 C is 
actually a maximum. In principle, the relevant "duality breaking" parameter (to 
be found by comparing the theory with the available data) is the discontinuity 
Ak' of the first derivative of /c e //(/3), namely 

Ak' = (V«(/3) - dpk%V3j) \^ • (21) 

However, from the practical point of view, it is much easier to work with the 
parameterization (|16| - (ll7p . where a non- vanishing Ak' is seen to be given by a 
non-vanishing Ad t = d\ + ^ — d\ ' ^ 0, i = 0, .., 3, which implies Ak' 0. 

Thus, some of the parameters d[ and d^ ' can be appropriately chosen in 
order to improve the agreement with the experiments. 

The parameterization p^|) and (jTTJ) is actually the simplest one being com- 
patible with the strong duality breaking, inspired by the exact solution of the 
two-dimensional model. So, in principle, the results obtained by this method 
could be improved if it were possible to explore the whole space of parameters 
more exhaustively. 

Because of limitations in computational resources, we are able to explore only 
in a very restrictive region of the space of parameters of [10) . Even in this case, 
the results are remarkably good if we make use of some hints. We will proceed 
as follows: first, we will chose some of the parameters to be equal to the best 
values found in [10] . Then, the novelty here will be to introduce a parameter 
that controls the strong duality breaking. In principle, one would expect the 
optimal set of parameters to be far from the set found in [10] . where quite 
simplifying hypothesis of minimal duality breaking were assumed. However, as a 
very encouraging signal, one finds that a very good agreement with observations 
is found by studying this region of the moduli space. Thus, from now on, we 
will fix 

Co = Co = 1-9389, Cx = Ci = ~, C 2 = a = 1-9205, 

and we will test the strong duality breaking in a small neighborhood varying 
A, C 2 and Ak'. 

Afterwards it will be useful to express k^A in terms of the low temperatures 



variable u as follows 

(d[- ] + d [ - ] u^ 



. 4-> - 4-^ 
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Now, let us move on and study the high and low temperature regimes. 



5 High and low temperatures 

The idea of this section is to find the optimal set of parameters (£3, A*, Afc'*) in 
Eqs. ([9]) and (|2Tj) that reproduce as close as possible the available Monte Carlo 
data at high temperature (see 32J ) . A hyper-cubic lattice has been chosen in the 
parameters space (every point in the lattice representing a possible set of high 
temperature parameters), then the free energy © will be evaluated at every 
point of the lattice. The optimal choice of parameters will be the one minimizes 
the following deviation function which, to some extent, represents the deviation 
between the ansatz and the Monte Carlo data, 

50 

X (C, v, A) = J2 \4 C d' X) (M +Io- F H A T IC (ft) 

i 

where 

&-A_i = ^ , [3 50 = (3 max = 0M, Jo = 2.4819 

f ht C (P) = 3cosh/3+ (3(tanh/3) 4 + 22(tanh/3) 6 + 187.5(tanh/3) 8 
+ 1980 (tanh/3) 10 + 24044 (tanh/3) 12 + 319170 (tanh/3) 14 - 

Here, we keep the terms up to the 15 th order of [32] since our algorithm is 
not sensitive to higher order terms. F H ^ lB the high temperature Monte 
Carlo free energy, (3 m can be assumed to be of ordeiQ 0.03, and I is a constant 
introduced for numerical convenience. 

In order to compare the Regge coefficients with the high temperature coef- 
ficients in , we need to change variable from to t = tanh ft 

Candidates to be the best parameters are 

C,; = 0.3273, = 0.1095, (23) 

while the discontinuity at the critical point is given by Afc' = 0.75 • 10~ 4 (which 
according to our precision is compatible with zero). As expected, the discontinu- 
ity of second derivative Afc" is of order one. This can be achieved, for instance, 
by the following choice of d^: 

4 +) = 1 , d[ +) = 1 , 4 +) = -0.1575 , 4 +) = 0-7116 
= 0.27 , d{~ } = 1 , = 0- 3498 > d V = 0- 6251 



7 Above P ?» 0.05 is not in the high temperature regime as the critical temperature is at 
j3* ss 0.22 which is only a factor of four larger. Indeed, (3 0.03 appears to be not small 
enough. Nevertheless, we will see that the agreement of our semi-analytical free energy with 
Monte Carlo data is excellent up to /3 ~ 0.03. 
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The agreement appears to be excellent: the deviation at high temperatures 
turns out to be <Jht(C*> v* , A*) ~ y/x(C > f*, A*)/50 » 3- 10~ 6 , which is certainly 
compatible with Monte Carlo results (see for instance pQ and references therein) . 

Once the optimal set of parameters have been found at high temperature, it 
is observed that such set also leads to a good agreement at low temperatures. 
The internal energy at low temperature is (k^J (/?) has to be expressed in terms 
of u as in Eq. ([22])) 

§) K \u) + 2I 1 =2u^F^ Ak '\u). (24) 

This is the average energy per spirjf] to be compared with (]f ) MC , the poly- 
nomial form in u which represents the Monte Carlo average energy for spin for 
small u found in [53] . 

We learned in [TU] that it is not convenient to use the expression on the 
right hand side of ([24)1 . Instead, the polynomial expression of [34] can be used 
to obtain the Monte Carlo estimation for the free energy at low temperatures; 
namely 

E\ KL _ d 



% ) («) + 2/i = 2u^-F LT MC {u) = £ af V 

' 8=6 



where 



F. 



LT 



MC 



Low temperature test function reads 

50 
i=l 

0.03 

Ui-Ui-i = — -, ui=u min , u 50 = M max = 0.03, 
50 



p MC 



(a) 



1 (12 6 60 10 
- — (it) H (it) 

2 V 6 v ' 10 w 

84 12 420 ,i4 

"i2 (u) + TT (u) ■ 



2i log u 



where Ji = —1 and J 2 = 0.0954. Our algorithm is sensitive up to the 15th order 
of the polynomial expression of [34], so we keep all these terms. M ma x has to be 
much smaller than u cr n = cxp(— 2[3 crit ) ~ 0.6, so it can be reasonably assumed 
to be of order 0.03. 

The deviation at low temperature between the Kallen-Lehmann and the 
Monte Carlo free energies evaluated for the same optimal parameters in (|23p . 



'to be more precise, both expressions can differ by a I\. 
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which have been found by asking the optimal agreement at high temperature, 
is CTLT (C*, v\ A*) « y/xLT(C, »*, A*)/50 « 7 • KT 6 . 

Remarkably, one finds the agreement at low temperature to be an order of 
magnitude better than that in the case of minimal duality breaking studied in 

Moreover, in the next section we will show that, besides improving the agree- 
ment with Monte Carlo results at low and high temperature, the implementation 
of the strong duality breaking also permits to describe features at the critical 
point. 



6 The Critical Point 

Once the parameters have been fixed as in (|23p one may verify that the behavior 
at the critical point is correctly reproduced as well. To do this one can fit near 
the critical point the non-analytiqf| part of the free energy in Eq. {§]) with the 
optimal parameters with a function of the form 

Fcrit « A + \f3~l3*\ 2 ~ a + c (3<f3 c 
F CRIT » A^\(3-(3*\ 2 - a +c (3>(3 C 

and find the optimal values of the constants c, A + , A_ and a, so that a will 
be our prediction for the critical exponent and A+/A- will be our estimate for 
the universal amplitude ratio. This form of the free energy's critical part is 
expected both from Conformal Field Theory and from experiments. The results 
of the fit doncP^I with Mathematica by fitting the (non-analytic part of the) 
Kallen-Lehman free energy with the optimal parameters in Eq. (I23|) with the 
above functions ^ from (3 mm = 0.22104 to (3 max = 0.22208 yield 



A_=4.18, A+=2.27, c = -0.19, a = 0.11 (25) 

the agreement appears to be very good when compared with recent estimations 
in [35], where the value a b s ~ 0.114(6) was found. In turn, we find 

Aa « a ° bs - a w 0.035 
a bs 

\ A +I A -)obs 

Both theoretical and experimental determinations of a can be found in [T]. 
The deviation Aa of our prediction (|25[) from observations appears to be less 
than 4% in all the more recent values, improving the minimal duality breaking 
results [10]. Figure ([lj shows the matching between the (non-analytic part of 



9 That is, one has to exclude the term log 2 cosh f) which does not contribute to the critical 
behavior. 

10 A regular sampling is chosen with step A/3 = 2* 10 -5 , thus the points of the sampling of 
Kallen-Lehmann are of the form P min + nA/3 and j3 max — mA/3 up to the critical point. 
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Figure 1: Bold points in the graph represent a sampling of Kallen-Lehmann free 
energy F 3 q ' (/3), while the continuous line corresponds to the critical part of the 
expected free energy Fcrit versus /3. One can observe the different leanings in the 
left and right sides of the critical point, which is ultimately related to a non-trivial A. 
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the) Kallen-Lehmann free energy and the Fcrit, and it confirms that in the 
parts of the graph of Fcrit in which the dependence on a is important the 
agreement is satisfactory. 

We expect that more general strong duality breaking can improve the very 
good agreement with observations and numerical data we found at high, low 
temperatures as well as at the critical point. For sake of simplicity (and be- 
cause of our computational resources) we considered the simplest form of strong 
duality breaking in which (inspired by the Onsager solution in two dimensions) 
the functional form of k e ff at high and at low temperature is the same (see 
Eqs. P^|) and (fT7]) h In turn, in our approach the cusp in the Fisher curve 
was implemented by the difference between the coefficients d[ + ^ at high 
temperatures and the corresponding coefficients d\ ^ at low temperature. How- 
ever, one could argue that it is likely the case the functional form of k e ff at 
low temperature in Eq. (|17p to be substantially different from the one at high 
temperature in Eq. (|16|) . In fact, there are infinite possible terms compatible 
with the Marchesini-Shrock symmetry that one can add to (fT7|) (still satisfy- 
ing conditions (fT8|) and (|20|) ) so that, in a sense, the very good agreement we 
found here is unexpected. It would be very interesting to explore theoretical 
arguments to further constrain k e ff. As a matter of fact, one may also try 
to add terms to k e ff compatible with the above physical requirements and see 
what happens. However, without a theoretical guide, the time required for an 
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exhaustive exploration of the moduli space would be out of range. The lack of 
sophisticated computational resources restricted us to explore a little region of 
the space of parameters. In particular, some parameters (Co, C2 > ^) were here 
kept fixed to take the same values found in [TO] just as a working hypothesis, 
and it is likely the case the optimal values to be slightly different from those 
of the ansatz in [TUj. Nevertheless, the implementation of the strong duality 
breaking led us to a quality leap, bringing our predictions at low temperature 
closer to Monte Carlo simulations. Remarkably, at the same time our model 
permits to reproduce the universal amplitude ratio A and to predict the critical 
exponent a. 

7 Conclusions 

In this paper we studied the distribution of Fisher zeros in the Kallen-Lehmann 
approach to the three-dimensional Ising model. It was shown that non-trivial 
amplitude ratios A = A + /A^ are compatible with the ansatz proposed in [3J. 
We proposed a mechanism (which we called strong duality breaking) by gener- 
ating a cusp in the curve of the zeros of the free energy F(f3) in the complex 
/3-plane, being the cusp located at the critical point. 

This mechanism not only permitted us to reproduce the Monte Carlo predic- 
tion for the value of A, but also led us to improve the results of [TU] bringing the 
predictions of the Kallen-Lehmann ansatz at low temperature closer to Monte 
Carlo simulations. The agreement turns out to be remarkable. In particular, 
the matching between our results and Monte Carlo estimations for the critical 
exponent a and for the amplitude ratio A exhibit a relative deviation of 3.5% 
and 7% respectively. 

We interpret our result as a motivation to continue the investigation of phe- 
nomenological semi-analytic methods of this kind as a promising approach to 
solve interesting problems in statistical physics. 
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